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Abstract We present a multiwavelength study of the 
formation of massive stellar clusters, their emergence 
from cocoons of gas and dust, and their feedback on 
surrounding matter. Using data that span from radio to 
optical wavelengths, including Spitzer and Hubble ACS 
observations, we examine the population of young star 
clusters in the central starburst region of the irregular 
Wolf-Rayet galaxy IC 4662. We model the radio-to- 
IR spectral energy distributions of embedded clusters 
to determine the properties of their HII regions and 
dust cocoons (sizes, masses, densities, temperatures), 
and use near-IR and optical data with mid-IR spec- 
troscopy to constrain the properties of the embedded 
clusters themselves (mass, age, extinction, excitation, 
abundance). The two massive star- formation regions 
in IC 4662 are excited by stellar populations with ages 
of ~ 4 Ma and masses of ~ 3 x 10 5 M (assuming a 
Kroupa IMF). They have high excitation and sub-solar 
abundances, and they may actually be comprised of 
several massive clusters rather than the single mono- 
lithic massive compact objects known as Super Star 
Clusters (SSCs). Mid-IR spectra reveal that these clus- 
ters have very high extinctions, Ay ~ 20 — 25 mag, and 
that the dust in IC 4662 is well-mixed with the emitting 
gas, not in a foreground screen. 



1 Introduction 

Massive stars in starbursts are thought to form pre- 
dominantly in Super Star Clusters (SSCs), which are 
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massive, compact, young stellar clusters that are candi- 
dates for present-day analogues of young globular clus- 
ters because of their similar sizes and inferred masses 
(e.g. [Whitmore et al.|[l993| |Schweizer et al.||1996[ ). Al- 
though the great majority of SSCs could be short-lived, 
destined to dissolve within ~ 10 megayears (Ma) and 



merge into the field star population (Fall et al. 2005 



Mengel et aL]|2QQ5| ) ? it is clear that their massive stars 



dominate feedback in starbursts. They return enriched 
matter and energy to the interstellar medium (ISM) via 
stellar winds and supernovae, and their hard radiation 
fields heat and excite gas and destroy molecules and 
dust grains. A high surface density of star formation 
produces galactic-scale superwinds that can expel the 
ISM from a galaxy (e.g. Heckman et al. 2000). 

SSCs are probably formed within giant molecular 
cloud complexes, so observations that can probe embed- 
ded dusty regions are required to detect the youngest 
SSCs. By the time SSCs have blown away enough 
of their natal material to be detected at visible wave- 
lengths, they are usually at least several million years 
old, so their most massive individual stars are expected 
to have already evolved off of the main sequence. At 
this stage it is difficult to infer much about SSC forma- 
tion environments, because the massive stars have sig- 
nificantly altered them. Infrared and radio wavelengths 
are best suited to studying the youngest SSCs, and ra- 
dio interferometry has revealed a population of heavily 
obscured, compact thermal radio sources with inferred 
nebular densities n e > 10 3 cm -3 and radii of a few par- 
sees that appear to be powered by 100s- 1000s of mas- 
sive stars. These sources were dubbed ultra-dense HI- 
IRs (UDHIIs, Kobulnicky and Johnson|1999[ ) in analogy 
with Galactic ultra-compact HIIRs (UCHIIs, e.g. Wood 



and Churchwell 1989), which are excited by one or a few 
massive stars on subparsec scales, with n e > 10 4 cm -3 . 

UDHIIs are thought to be young, embedded SSCs 
that have not yet disrupted and expelled the majority 
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IC4662 
D = 2Mpc . 
10pc~l M 




IRAC 3.3, 5.8, 8 pm 
ATCA 3cm contours 
(Johnson at al. 2003) ■ 




Fig. 1 IRAC image of IC 4662 (with red, green, blue in 
3.3, 5.8, 8 /xm bands respectively) features three dusty re- 
gions that are coincident with the optically identified star- 
formation regions A2, Al, and B (from left to right) of 



Heydari-Malayeri et al. ( 1990 ). Bright thermal radio sources 



are shown in ATCA 3 cm contours (levels start at 3cr, steps 
are 3cr, |Johnson et al.|2003| 




of their natal and surrounding material. Their hard 
intrinsic emission is absorbed and reprocessed by sur- 
rounding dense gas and dust to emerge as a thermal 
blackbody component in the mid-IR, where it peaks 
at ~ 60 — 100 fim depending on the temperature and 
distribution of the dust (e.g. Vacca et al. 2002). In 
the radio regime, UDHIIs are identified by a thermal 
bremsstrahlung spectrum that is self-absorbed below a 
few cm (e.g. | Johnson et aL||2003 ). These features dom- 
inate the spectral energy distribution (SED) of UDHIIs 
and can be used along with any escaping visible light 
emission to constrain the temperature, density, and ge- 
ometry of the gas and dust excited by the SSC, as well 
as the age and stellar content of the cluster itself (e.g. 
Martin-Hernandez et al. 2005 2006). Spectroscopic 



measurements provide further constraints on these pa- 
rameters as well as the excitation and composition of 



the emitting material (e.g. Verma et al.||2003 ). 



2 Spitzer Imaging of Embedded Massive 
Star-Formation Regions in IC 4662 

IC 4662 is a nearby (d « 2 Mpc, 1" ~ 10 pc), iso- 
lated dwarf galaxy that contains two bright Wolf-Rayet 
(WR) star-forming regions that are extended over ~ 50 



pc, and very little molecular gas and dust ( Heydari^] 



Malayeri et al.|1990 ). Radio observations at 2 and 6 cm 
revealed that compact thermal sources containing up to 
a few hundred massive stars are associated with the WR 




Fig. 2 Top: ACS image of IC 4662 (red, green, blue in 
F435W, F658N, F814W, respectively) shows regions A2 and 
Al overlaid with 3 cm radio contours. Ha and radio con- 
tinuum emission have similar overall morphologies, but Ha 
has filaments, sharp edges and apparently empty regions. 
Bottom: Overlay of 8 /xm contours shows that warm dust 
emission peaks where Ha peaks 



star-forming regions, but they are not centered on the 



brightest optical clusters within them (Johnson et al. 



|2003p . We obtained Spitzer IRAC and MIPS imaging 
and IRS high-resolution spectroscopy of these regions 
in order to explore the mid-IR properties of the stars, 
gas, and dust in these regions. We also obtained Hub- 
ble Space Telescope (HST) ACS imaging to study the 
optical sources in these regions at high resolution, but 
we focus on Spitzer data in this contribution and treat 
the starburst regions as unresolved. 

A three-color image of IC 4662, constructed from 
the IRAC 3.3, 5.8, and 8-//m bands, is shown in Fig- 
ure [I] along with an overlay of the 3 cm radio contours 
from | Johnson et ah (2003). The galaxy has three bright 



extended regions in the IRAC bands, two of which cor- 
respond to the radio and optical regions Al and A2, 
and a third (region B, not included in the HST field 
of view) that is offset to the southwest and fainter in 



SSCs in IC 4662 
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Fig. 3 IC 4662 regions Al (black, blue) and A2 (red) 
SEDs and fits to a model of an SSC embedded in a UDHII 
embedded in a dust cocoon, plus a stellar SED at A < 5 /mi 

all bands. Figure [2] compares both the radio and IRAC 
8 /im contours with the optical emission on an ACS 
three-color image (in filters F435W, F658N, F814W). 
The radio contours follow the extended Ha morphol- 
ogy, and the extended 8 /im emission, which is a tracer 
of PAHs or hot dust depending upon metallicity (e.g. 
Engelbracht et al.||2005| |Wu et al.|[2QQ6] ) ? peaks at the 
Ha peaks. Regions Al and A2 are extended on scales 
of ~ 50 pc = 5" , but they are complex and may contain 
sources that are too compact to be resolved by Spitzer's 
smallest pixels, which are 1.2" for IRAC. 

We performed photometry on regions Al and A2 
in the Spitzer pipeline-reduced IRAC and MIPS data. 
The apertures were matched across all Spitzer bands at 
a radius of ~ 7" , for which we also recomputed ATC A 
radio fluxes for consistency. The IR background was 
subtracted using Spitzer SPOT predictions. In order 
to incorporate the lower-resolution IRAS measurements 



(Beichman et al. 1988b into our analysis we scaled the 



IRAS fluxes so that the 25 /im flux matches a power- 
law fit to the Spitzer 24 and 70 /im data interpolated 
at 25 /im. This yields three additional independent es- 
timates of fluxes of Al and A2 at 12, 60, and 100 /im. 
The resulting SEDs extend from 3 cm to 3 /im and fea- 
ture a radio free- free component longward of about 150 
/im, a warm dust blackbody peaking at around 100 /im, 
and the onset of stellar photospheric emission appears 
in the 3 /im point (Figure [3]). 




Fig. 4 Constraints on age r and mass M for the exciting 
clusters in IC 4662 regions Al and A2: radio flux requires 
ionizing photon rate Q > 10 51-4 s _1 (above green line) and 
bolometric luminosity Lboi > 10 7 L© (above blue line); lack 
of observed 3 cm non-thermal emission implies that the clus- 
ters are outside of the SNe box (to left / right of or below red 
line); fit to the stellar part of the SED (A < 5/im) requires M 
~ 3 — 30 x 10 5 M (at ages on the magenta line); detection 
of optical and WR lines implies ages of 3 — 5 Myr (between 
the yellow lines). In summary, the exciting clusters of Al 
and A2 must have ages of about 4 Ma and masses of about 
3 x 10 5 M 



3 Modeling the SEDs of UDHIIs 

In order to interpret the measured SEDs, we begin by 
considering a simple modeQ of a star cluster embed- 
ded in an HII region which is embedded in a dust co- 

fl2002D and motiv a ted by 



coon, following 



Vacca et al. 



Churchwell et al. (1990) 



the UCHII models of e.g. 
The SED produced by this system is the sum of an op- 
tically thick, radio free-free component (assumed to be 
a constant-density, dust-free HIIR at a temperature of 
10 4 K extending from the center to an inner radius) and 
an optically thin dust shell component (extending from 
the inner to an outer radius with a constant density 
of silicate grains and a power-law temperature profile). 
We fit the SEDs from 3 cm to 5 /im of Al (A2) to this 
model and arrive at HIIR parameters: half-light radii 
of 4 (2.4) pc, n e = 1.2 (2.6) x 10 3 cm -3 , and gas masses 
of about 2 x 10 5 M . The dust shells have tempera- 
tures at the inner radius of 465 — 485 K and masses 
of about 10 3 M . The low model gas and dust masses 



relative to other UDHIIs like those in He 2 — 10 (Vacca 



1 We will consider a more complete and complex model such as 
DUSTY jNenkova et al. |2000|> in future work, as we conclude in 
§ |4.1| that IC 4662 requires a mixed model for its dust rather than 
a foreground screen. 
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Fig. 5 IRS slit nod positions overlaid on IC 4662 regions 
Al and A2, shown on ACS F658N (left) and IRAC 8 fim 
(right). The smaller SH slits are perpendicular to the larger 
LH ones 



et al. 2002 ) suggest that the exciting clusters Al and A2 
are more evolved, and have destroyed or ejected more 
of their surrounding matter in the process of emerging 
from their cocoons. 

Figure [3] shows that the HIIR+dust cocoon model 
fits the data points well, with the exception of the near- 
IR filters that are dominated by stellar continuum emis- 
sion and may contain PAH emission at 3.3/im. The 
radio data provide several constraints on the age and 
mass of the ionizing stellar clusters within regions Al 
and A2 that are illustrated in Figure [4] From the radio 
continuum fluxes we derive a Lyman continuum pho- 
ton rate Q ~ 2 — 3 x 10 51 s -1 , which requires at least 
a certain mass in massive stars and hence an overall 
cluster mass for a full initial mass function as a func- 



tion of cluster age r (we adopt a Kroupa IMF, Kroupa 



2001). Clusters in the M— r space above the green line 
in Figure [4] have enough massive stars to explain the 
observed Q values even if the radio fluxes are under- 
estimates due to self-absorption. The blue line shows 
a weaker constraint following the same argument for 
the bolometric luminosity Lboi of the exciting cluster. 
The lack of non-thermal emission at 3 cm implies no 
observable supernova (SN) activity, which excludes the 
region in M— r that is inside the red box. The detection 



of He II 468.6 nm line emission ( Heydari-Malayeri et al. 



1990) further restricts the permitted age range for the 
exciting sources to lie between 3 and 5 Ma, as shown 
by the yellow lines. Finally, a fit to the 3.3 and 4.5 /xm 



points in the SEDs of Bruzual and Chariot (2003) mod 



els for a range of cluster ages and masses yields r and 
M values that fall on the magenta points, so the allowed 
cluster parameters are very tightly constrained at r ~ 4 
Ma and M « 3 x 10 5 M . Thus we have added a stel- 
lar component to our SED model from the appropriate 
Bruzual and Charlot| ( |2003| ) model in Figure |tj 
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Fig. 6 IRS SH spectra of IC 4662 regions Al and A2 
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Fig. 7 IRS LH spectra of IC 4662 regions Al and A2 

4 Extinction, Excitation, and Abundances 
from Spitzer IRS Spectroscopy 

We observed IC 4662 regions Al and A2 with Spitzer 
IRS in both high-resolution (~ 500 km s _1 ) settings: 
short-high (SH) covers about 10 — 20 fim with a ~ 

5 x 11" slit, and long-high (LH) covers 19 — 37 /im with 
an ~ 11 x 22" slit. Both slits are shown in both nod 
positions for each object in Figure [5] on the continuum- 
subtracted Ha and 8 /im images. The spectra were 
reduced and extracted from the full apertures using 
SPICE, and background-subtracted using the SPOT 
background calculator. Figures [6] and [7] show the SH 
and LH spectra, which feature a warm dust continuum, 
many high-excitation fine-structure lines, a weak PAH 
feature at 11.3 /im but no discernible silicate absorp- 
tion, and broad H I Humphreys a (Hua) emission at 
12.37 /im. Fine-structure lines in the IRS spectra make 



SSCs in IC 4662 
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possible the determination of many properties of the 
gas, dust, and stars within the emitting regions. 

4.1 Extinction and Electron Density 

We can estimate the extinction by comparing theoreti- 
cal ratios of fine-structure lines that are well-separated 
in A with the observed values. In this case we used 
[Ne ill] 15.6/36 /im, which does not vary substantially 
over a wide range of expected electron temperatures 
and densities (T e , n e ). We thus assume T e = 10 4 K and 
n e = 10 2 cm -3 and use the mid-IR extinction law for 



the Small Magellanic Cloud of Weinga rtner and Draine| 
( |2QQ1[ ) to derive for Al (A2) an extinction A v = 11 (9) 
mag for a foreground dust screen or Ay = 25 (20) mag 
for dust mixed with the emitting gas. For either dust 
geometry, correcting the observed [S ill] 18.7/33.5 jam 
ratio for this extinction leads to n e « 150 (212) cm -3 
for Al (A2), so our initial assumption of a low n e gives 
a consistent result, whereas a higher assumption gives 
an inconsistent on^| Our fit to the SED (§ [5]) yields a 
higher n e ~ 10 3 cm -3 , which might be due to aperture 
differences between SH and LH slits: the larger area 
viewed by the LH slit could produce larger fluxes for 
the LH lines and hence drive down the measured ratio 
and derived n e . We also note that Hunter et al. (2001) 



reported a high gas density in warm photo-dissociation 
regions of n = 2.8 x 10 4 cm -3 from ISO data. These 
differing density tracers probe different gas phases that 
may not share the same emitting volume. Here we 
adopt a low value of n e = 10 2 cm -3 in order to ob- 
tain a lower limit on the extinction. 

These values for Ay greatly exceed those that were 
determined based on optical spectroscopy by |Heydari-| 
Malayeri et all (119901), who reported A v = 0.27 



mag from B aimer decrement measurements, or even 
Stasihska et al. (1986), who reported a much higher 
optical extinction equivalent to Ay = 2.5 mag. Both 
groups used apertures of 4—5" FWHM, which is compa- 
rable with the spatial extent of the SH aperture, which 
encompasses the bulk of the mid-IR and Ha emission 
from Al and A2 (Figure [5]), so comparison of our data 
sets is warranted. Only a mixed dust model can rec- 
oncile the low optical Ay values with the high mid-IR 
ones. If the extinction were as low as that inferred from 
visible wavelengths, the mid-IR line ratios would have 
the theoretical values. If the dust were in a foreground 
screen, the high IR extinction would obscure most of 
the optical recombination line emission; however, we 

2 By comparing ratios for [Ne ill] and [S ill] we are assuming 
they are co-located, although their ionization potentials do differ 
significantly. 



detect only about three times the optically-derived Ly- 
man continuum flux from mid-IR and radio data, so a 
screen geometry is ruled out. In other Spitzer obser- 
vations of nearby and hence well-resolved massive star- 
formation regions, dust has been observed to be mixed 



with the stars and gas (e.g. RCW 49, Churchwell et al. 



2004). 



4.2 Excitation and Massive Star Census 

Several of the mid-IR fine-structure line ratios that we 
can measure from the IRS spectra are useful excitation 
indicators for star- forming regions. For example, star- 
bursts and HIIRs form a diagonal band in a plot of [S iv] 
10.5/[S ill] 18.7 /im versus [Ne ill] 15.5/[Ne n] 12.8 /im, 
where higher-excitation objects have higher values of 
both ratios (e.g. |Wuet al.|2006[ ). With [S iv] 10.5/ [S III] 
18.7 /im ratios > 1 and [Ne ill] 15.5/[Ne n] 12.8 /im ra- 
tios > 7, IC 4662 Al and A2 both fall in the upper 
right corner of such a diagram, along with other metal- 
poor WR systems like NGC 5253 and II Zw 40 (e.g. 
Verma et aLl|2003| ). |Heydari- Malayeri et al] fll990| ) also 
report a very high excitation for IC 4662 based on its 



[Ne ill] 386.9 nm/H/3 ratio, and Mas-Hesse and Kunth 
(1999) report a high effective temperature T e ff ~ 38000 
K from its UV spectrum. 

High excitation is caused by the hard radiation field 
of the most massive stars that are present in SSCs, and 
we have several handles on the massive star content of 
IC 4662's star- forming regions from our data. Assuming 
that the radio continuum flux is due to free-free radia- 
tion from HIIRs as above, we find Q = 4.6(3.0) x 10 51 
s _1 for Al (A2), and from the Hua flux we find Q 
= 3.1(2.6) x 10 51 s _1 , in reasonable agreement. From 
these Lyman continuum luminosities we infer that the 
number of hot stars needed to power the HIIRs is 
- 260 - 460 07V stars. From e.g. |Heydari-Malayeri| 



et al. (1990) it is known that these regions contain WR 



stars, so using the models of Schaerer and Vacca (1998) 
and Starburst99 QLeitherer et al.| 19*99 ) we estimate that 
~ 70 — 150 WRs should be present. This number agrees 



with that of Hey dari- Malayeri et al. (1990) after cor- 



recting for the mid-IR-derived extinction. 
4.3 Abundances 

We calculated neon and sulfur abundances for regions 
Al and A2 from the observed fluxes of fine-structure 
lines relative to that in Una at 12.37 /im, following a 



procedure similar to that used by Wu et al. (2007) 



whose IRS abundance study of blue compact dwarf 
galaxies (BCDs) is useful to compare with ours. We 
measured sulfur in two excited states at [S ill] 18.71 
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fim and [S iv] 10.51 /im, which we expect to comprise 
most of the sulfur at such high excitation. Following 
Wu et al.| ( |2QQ7| we assumed that only 10 % of the 
total S is in [S II], since the typical ionic abundance 
of [S ii] in BCDs is 10 - 15 % of that of [S ill] ^ 
1997| ). For neon we used the fluxes of 



tov et al. 1994 



[Ne n] 12.81 /im and [Ne ill] 15.55 /im and made no cor- 
rection for unseen ionization states. All of these lines 
were observed with the same slit (SH) as Hua so that 
no aperture corrections are needed, as would be needed 
if we used the LH detections of [S ill] 33.48 /im and 
[Ne ill] 36 /im. 

The observed fluxes were corrected for extinction us- 
ing the mixed model results discussed in § |4.2[ and 
the resulting abundances for Al (A2) are: Ne/H = 
1.4 x 10~ 5 (2.0 x 10~ 5 ) and S/H = 2 x 10~ 6 (2.4 x lO" 6 ). 
In order to compare directly with the BCDs of |Wu et ah] 
( |2007 ), we adopt their choices for solar abundances of 
(Ne/H) = 1.2 x 10- 4 ( iFeldman and Widing|2003l ) and 



(S/H) = 1.4 x 10~ 5 ( |Asplund et al.||2005[ ). Thus in 
solar units we have (Ne/H) = 0.12(0.17) (Ne/H) and 
(S/H) = 0.14(0.17) (S/H) . |Heydari-Malayeri et al.| 
(1990) report a mean Ne/H abundance of 2.25 x 10 -5 , 
which is consistent with our lower limits but nearly a 
factor of two larger. We note that the Una is signifi- 
cantly broader than the other nebular lines in the IRS 
spectra. We investigated whether this broadening could 
be the direct signature of WR stars via their He i/ll 
12.36 /im emission ( |Smith and Houck||2001[ ). However, 
the number of WR stars needed to explain the puta- 
tive excess Una flux is far larger than the number of 
O stars expected to be present. Therefore the cause of 
the broad Hua line remains uncertain. 

The relatively weak PAH feature at 11.3 /im is con- 
sistent with the trends found by |Wu et aL ( 2006 ) for 
a sample of BCDs that have high excitation and low 
abundance like IC 4662 regions Al and A2. 



5 Conclusions 

The nearby dwarf irregular galaxy IC 4662 harbors two 
sites of recent massive star formation that have thermal 
radio spectra, bright dust and Ha emission, and rich 
mid-IR spectra that feature nebular fine- structure lines, 
a hint of PAHs at 11.3 /im, and broad Hua emission. 
Our emission-line analysis indicates that this starburst 
has high excitation and low abundances like many UCD 
galaxies and in agreement with literature values. How- 
ever, we find much higher extinctions (Ay ~ 20 — 25 
mag) than previous shorter- wavelength studies, which 
can only be reconciled with the optical observations by 
a mixed geometry for the gas and dust in these re- 
gions. SED fitting of the radio-to-near-IR SEDs of the 



UDHIIs in regions Al and A2 suggests that they are 
more evolved toward the SSC stage than similar em- 
bedded objects (e.g. in He 2 - 10 and NGC 5253): Al 
and A2 have lower gas masses in their HIIRs and dust 
masses surrounding them, and this is consistent with 
their older ages of about 4 Ma, inferred from radio and 
optical/near-IR data. The clusters that power the UD- 
HIIs have masses of about 3 x 10 5 M (assuming a full 
Kroupa IMF), but they may not be massive, compact, 
monolithic SSCs like those found in many starbursts, 
but rather clusters of 10 3 — 10 4 M clusters that fill a 
larger volume of ~ 50 pc diameter. 
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